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THE ARCTIC CIRCLE BOUNDARY AND THE AIRY PROCESS 

By Kurt Johansson 1 

Royal Institute of Technology 

We prove that the, appropriately rescaled, boundary of the north 
polar region in the Aztec diamond converges to the Airy process. 
The proof uses certain determinantal point processes given by the 
extended Krawtchouk kernel. We also prove a version of Propp's con- 
jecture concerning the structure of the tiling at the center of the Aztec 
diamond. 

1. Introduction and results. Domino tilings of the Aztec diamond were 
introduced in [8, 9]. Asymptotic properties of random domino tilings of the 
Aztec diamond have been studied in [5, 12, 15]. In particular, in [12] the 
existence of the so-called arctic circle was proved. The arctic circle is the 
asymptotic boundary of the disordered so-called temperate region of the 
tiling. Outside this boundary the tiling forms a completely regular brick 
wall pattern. The methods in [12] combined with the results in [13] show 
that the fluctuations of the point of intersection of the boundary of the 
temperate region with a line converge to the Tracy-Widom distribution of 
random matrix theory. In this paper we extend this result to show that the 
fluctuations of the boundary around the arctic circle converges to the Airy 
process introduced in [23]. The paper is a continuation of the approach used 
in [14] and [15], where certain point processes with determinantal correla- 
tion functions [24] and the Krawtchouk ensemble, were used. We will use 
the general techniques developed in [16] and investigate an extended point 
process which also has determinantal correlation functions given by a kernel, 
which we call the extended Krawtchouk kernel. 

The Aztec diamond, A n , of order n is the union of all lattice squares 
[m, m + 1] x [1, 1 + 1], m, I € Z, that lie inside the region {(x%, y%); \xi\ + \yi\ < 
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n + 1}. A domino is a closed 1 x 2 or 2 x 1 rectangle in M 2 with corners in 
Z 2 , and a tiling of a region C M 2 by dominoes is a set of dominoes whose 
interiors are disjoint and whose union is R. Let T(A n ) denote the set of all 
domino tilings of the Aztec diamond. The basic coordinate system used here 
will be referred to as coordinate system I (CS-I). 

Color the Aztec diamond in a checkerboard fashion so that the leftmost 
square in each row in the top half is white. A horizontal domino is north- 
going (N) if its leftmost square is white, otherwise it is south-going (S). 
Similarly a vertical domino is west-going (W) if its upper square is white, 
otherwise it is east-going (E). Two dominoes are adjacent if they share an 
edge, and a domino is adjacent to the boundary if it shares an edge with the 
boundary of the Aztec diamond. The north polar region (NPR) is defined 
to be the union of those north-going dominoes that are connected to the 
boundary by a sequence of adjacent north-going dominoes. The south, west 
and east polar regions are defined analogously. In this way a domino tiling 
partitions the Aztec diamond into four polar regions, where we have a regular 
brick wall pattern, and a fifth central region, the temperate region, where 
the tiling pattern is irregular. 

Let T E T(A n ) be a tiling of the Aztec diamond and let v(T) denote the 
number of vertical dominoes in T. We can define a probability measure on 
T(A n ) by letting vertical dominos have weight a and horizontal dominos 
weight 1, that is, 

a v(T) 

(i-i) m = m . 

In this paper we will prove the asymptotic results for the uniform case (a = 1) 
only to keep the asymptotic analysis simpler. 

We will study the part of the boundary of the NPR which lies above a 
neighborhood of x\ = 0. To define the boundary we will use certain nonin- 
tersecting paths which describe the domino configuration and which are also 
essential in our analysis below. They were called DR-paths ([25], page 277) 
in [15]. On a W-domino placed so that it has corners at (0,0) and (1,2), we 
draw a line from (0,1/2) to (1,3/2). On an E-domino placed in the same 
position, we draw a line from (0,3/2) to (1, 1/2), and on an S-domino placed 
so that it has corners at (0,0) and (2,1), we draw a line from (0,1/2) to 
(2,1/2). We do not draw any line on an N-domino. As discussed in [15], 
these lines will form nonintersecting paths from A r = (— n — 1 + r, — r + 1/2) 
to B r = {n + 1 — r, — r + 1/2), 1 < r < n. The top curve, from A\ to B\, 
can be viewed as a function, t — > X n (t), \t\ < n, in CS-I. We will call X n {t) 
the NPR-boundary process, see Figure 1. The NPR is exactly the part of 
the domino tiling that lies above X n (t), and consits only of N-dominoes. 
Between the nonintersecting paths there are other regions of N-dominos. 
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Before we formulate the limit theorem for X n (t), we recall the definition 
of the Airy process. The extended Airy kernel is defined by 

A(t,£;t',0 

{/•oo 
^ e- x ^- T )Ai(e + A)Ai(e' + A)dA, if r>r', 

-j e~ A ( r ~ r ') Ai(e + A)Ai(e' + A)dA, if t<t', 

for r, t' , £, £' £ E. Fix n < ■ • ■ < r m and let A m = {ri, . . . , r, n }. Let \i be the 
product of counting measure on A n and Lebesgue measure on M. Define / on 
A m x K by f(rj,x) = l(£j,oo) ( x ) f° r given numbers £1, . . . , £ m . It is proved in 
[16] that f 1 / 2 (T,x)A(T,x;a,y)f 1 / 2 (a,y) is the integral kernel of a trace class 
operator on L 2 (A m x M,/x). The Airy process, t^A(t), is the stationary 
stochastic process whose finite-dimensional distributions are given by 

(1.3) PL4(n) < £i, . . .,A(r m ) < £ m ] = det(J - / 1/2 ^/ 1/2 )^(A m xiR)- 

It is proved in [23], see also [16], that A(t) has continuous paths. The dis- 
tribution of A(t) is F2, the Tracy- Widom distribution for the largest eigen- 
value of a GUE-matrix [27]. It has recently been shown that the distribution 
function in (1.3) satisfies certain differential equations [1, 28]. 
Our main result is the following: 

Theorem 1.1. Let X n (t) be the NPR-boundary process and A(t) the 
Airy process as defined above and let the weight a in (1.1) be equal to 1. 
Then, 

as n — > 00, in the sense of convergence of finite- dimensional distributions. 

The theorem could be extended to a general weight a in (1.1) by the same 
type of argument. We restrict to a = 1 for simplicity. Similarly, it is possible 
to show convergence along other parts of the boundary, except right near the 
point where the arctic circle is tangent to the asymptotic square containing 
the tiling where the boundary behavior is different. See the remark in the 
last paragraph of Section 2. 

As was done in [16] for the convergence of the interface of a polynuclear 
growth model to the Airy process, it should be possible to extend Theorem 
1.1 to a functional limit theorem. It would then follow that maxtA^i), 
suitably rescaled, asymptotically has F\ fluctuations, where F\ is the Tracy- 
Widom distribution for the largest eigenvalue of a GOE-matrix [27]. 

The above result also has an interpretation as a convergence theorem for 
a certain polynuclear growth model, see [15], Section 2.4. The polynuclear 
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Fig. 1. An N PR-boundary process. 



growth model studied in [23] is a limiting version of this model (a — > 0, 
n — > oo at appropriate tates). The NPR-boundary is also directly related to 
the shape in the corner growth model studied in [13] and [15], Section 2.4, 
and, hence, also to the totally asymmetric exclusion process. If G(M, N) is 
the last passage time as in [13], the results of the present paper show that 
the fluctuations of the boundary of the shape fit = {(M,N) ; G(M,N) + 
M + N — 1 < t} close to the diagonal (N, N) , when q = 1/2, converges to the 
Airy process. This can be extended to the other parts of the boundary, away 
from the axes, and to other values of q. This extends the relation between 
the Meixner ensemble and the Krawtchouk ensemble given in Lemma 2.9 in 
[15], see also [20]. 

The theorem will be proved using a certain determinantal point process. 
One way of seeing this point process is to put a green particle at the center 
of the black square on each S-domino and W-domino, and a red particle at 
the center of the white square of each S- and W-domino. These dots define 
a point process and we will see below that it has determinantal correlation 
functions. It is directly related to the nonintersecting paths defined above 
and also to the zig-zag paths around black and white squares discussed in 
[8, 9, 14]. The precise definition of this point process will be given in Section 
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2. The asymptotic results needed to prove Theorem 1.1 will be discussed in 
Section 3. It has been conjectured in [12] that at the center of a uniform 
random tiling of the Aztec diamond the tiling looks like a random domino 
tiling of the plane under the Burton-Pemantle measure [4]. A proof of a 
version of this conjecture will be given in Section 4. 

It is natural to conjecture that boundaries between irregular and regular 
tiling regions in many other two-dimensional random tiling problems are 
also described by the Airy process. In terms of the so-called height function, 
which we will not define here, this would mean that we would see the Airy 
process where we have a boundary between a flat surface and a curved 
surface. This type of result can also be proved in the rhombus tiling problem 
discussed in [22]. This has been done recently in [10]. In this problem it 
would also be possible to use the formulas derived in [22], which would 
lead to computations very similar to those in the present paper. Another 
candidate where it may be possible to prove convergence to the Airy process 
would be for fluctuations around the arctic circle (ellipse) in rhombus tilings 
of hexagons [7, 15], but here the asymptotic analysis appears to be more 
difficult [2]. 

2. The point process. Introduce a new coordinate system (CS-II) with 
origin at (— n, — 1/2) and axes en = (1,1), fn = ( — 1, 1), which gives the 
following coordinate transformation between CS-I and CS-II, 

(9U x 1 = x 2 -y 2 -n, 

{ ' yi = x 2 +y 2 -l/2. 

In CS-II the nonintersecting paths defined in Section 1 go from Aj = (0, — j + 
1) to Bj = (n + 1 — j, — n), 1 < j < n, and have three types of steps (1,0), 
(0,-1) and (1,-1), see Figure 2. These nonintersecting paths specify the 
tiling uniquely. The measure (1.1) is obtained by letting the steps (1,0) and 
(0,-1), which correspond to vertical dominos, have weight a and the step 
(1, —1) weight 1. 

To formalize this, let Q = (V, E) be a directed graph with vertex set V = 
N x Z and directed edges from to (i + — 1) and (i + — 

1) for i > 1, and from (0,j) to (1, j — 1) and (l,j). The edges from 
to (i + l,j — 1) have weight 1, whereas all other edges have weight a. A 
path 7r from A&VtoBGVisa directed path along succesive directed 
edges starting at A and ending at B. The weight uj{tt) of a path tt is the 
product of the weights of all the edges in the path. Two paths tt and tt' are 
nonintersecting if they do not share a common vertex. If A = (A\, . . . , A m ) 
and B_= (B\, . . . ,B m ), where A r ,B r E V, then V n .i {A;B_) denotes the set 
of all nonintersecting paths tt\ , . . . , 7r m , where ttj goes from Aj to Bj . The 
weight of a family 7r = (in, . . . , ir m ) of paths is uj(n) = YlJLi Take 
N > n and set Aj = (0, 1 — j) and Cj = (n, — n + 1 — j), 1 < j < N; see 
Figure 3. 
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Lemma 2.1. If (tt, . . . , 7T/v) are paths inVrL±(A;C_), thenirk goes through 
the vertex B\. = (n + 1 — k, —n), 1 < k < n. 

Proof. Since tt\ ends at C\ = B±, the claim is true for k = l. Each 7Tfc, 
1 < k < n, has to go through one of the points B — 2, . . . , B n since C2, ■ ■ ■ ,C n 




Fig. 2. CS-11 and nonintersecting paths describing the tiling. 




Fig. 3. The nonintersecting paths in the graph Q . 
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lie below the line yi = —n and all the paths are right/down directed. The 
nonintersection constraint implies that ir k +\ has to pass this line to the left 
of 7Tfc and the claim follows. □ 

By this lemma there are well-defined projections 

P:n = V n .iXA 1 ,...,A N ;C 1 ,...,C N )^n' = Vnx(A 1 ,...,A n ;B 1 ,...,B n ) 
and 

Q : Q — > $7" = Vn.i.(Bi, ■ ■ ■ , B n , A n+ i, . . . , A^; C±, . . . , CV), 

such that lj(tt) = uj(P(tt))uj(Q(it)). An event D in the domino tiling of the 
size n Aztec diamond corresponds, via the bijection, to an event D C Q' . Set 
D* N = {vr G 0; P(vr) G D}. Then, 



(2.2) 



E^n^U) 

(E 7r , e5 ^(E / ))(E 7r " e n^(2L / )) 

P[L>]. 



(E^n^U'OXE^w^U'O) 
The right-hand side in (2.2) is independent of N > n and, hence, 
(2.3) F[D} = hm P^ra. 

AT— >oo 

We want to map, bijectively and preserving weights, the nonintersecting 
paths in H to a new family of nonintersecting paths which is appropriate 
for the application of the general results in [16]. The idea is as follows. Each 
path Tij, from to has a first and a last vertex on each vertical line 
X2=j- The first and the last vertex could, of course, coincide. We want to 
put the first and the last vertex on different vertical lines. These first and 
last vertices form the point process we are interested in. Let Q' = (V',E') 
be the directed graph with V' = N x Z and directed edges from (2i — 1, j) to 
(2i,j), i > 1, from (2i,j) to (2i,j - 1), from (2i,j) to (2i + 1, j), horn (2i,j) 
to (2i + l, j + 1), i > 1, from (0,j) to (1, j) and from (0, j) to (l,j + l), j G Z. 
We put the weight a on the edges from (2i,j) to (2i + 1, j + 1) and from 
(2i,j) to (2i, j — 1), i > 1, J G Z. All other edges have weight = 1. 

We can describe a path tt^ from to by giving the first, Pk(j), and 
the last, Qk(j), point on each verical line x% = j, < j < n; Pfe(O) = Qfe(O) 
always. From Qk(j) *° -ffc(j + 1) we take either a step (1, 0) or a step (1, — 1), 
and from Pk(j) to Qk(j) we take a certain number, > 0, of down steps 
(0, -1). Map Q k (j) = (j,q) to R k (2j) = (2j,q + j), < j < n, and P k (j) = 
(j,p) to R k (2j-1) = (2j-l,p+j). A step from Q k (j) to P k {j + l) is mapped 
to two steps. One step from R k (2j) to R k (2j + 1), which is either (1,0) or 
(1, 1), and then a fixed step from R k (2j + 1) = (2j + l,p + j) to (2j + 2,p + j), 
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Fig. 4. TVie nonintersecting paths in the graph Q' corresponding to the tiling in Figure 
2. The particles in the determinantal process are the circled dots. 



that is, a step (1,0). The vertical steps from Pk(j) to Qk(j), if there are 
any, are mapped to the same number of vertical steps from (2j,p + j) to 
(2j,q + j). Set A' k = (2n, 1 — k), 1 < k < N. The above procedure maps 
V n .i. (A 1 ,...,A N ;Ci,..., C N ) to V n .i. (A ± , . . . ,A N ;A[, . . . ,A' N ) bijectively, see 
Figure 4. Also, the weight of a family of paths is preserved. The probability 
of a family of paths 7r in £1^ = Vn.i. (A]A!_) is 

(2.4) M.)- ^ 



N 



An event -Djy in as above maps to an event Dn in fi^r using the bijection, 
and from (2.3), we have F[D] = limjv-»ooPjv(^iv)' 

The paths 7ri, . . . , ttn are uniquely specified by the points (Rk{j))i<j<2n~i,i<k<N ■ 
We will write Rk(j) = (j,xQ, 1 < J < 2n — 1, l<k< N and x = (x 3 k ); x k = 
x| n = 1 — k, k > I; x\> X2> • • • > x~L. These points define a particle configu- 
ration in {1, . . . , 2n — 1} xZ and we obtain a point process on this space which 
we can think of as 2n — 1 copies of Z. We see that the NPR-boundary pro- 
cess is obtained by joining Qi(0), -Pi(O), -Pi(l), Qi(l), ■ • ■ , Pi(n), Qi(n) with 
straight lines in this order. If Q\{j) has coordinates (j,qj) in CS-II, then 
i?i(2j) = (2j,qj + j), so qj = x 2 ^ — j. Similarly, if P\(j) has coordinates 
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(2.5) 



(j,Pj) in CS-II, then pj = x 




j. We see from (2.1) that 



2j-i 



1/2), 



1/2) 



in CS-I. 



The red points described in the introduction are Pk(j) + (— 1/2, 0) in CS-I 
and the green points are Qk(j) + (1/2,0); note that we can have Qk(j) = 
Pk(j)- This is seen from how we defined the nonintersecting paths and the 
red and green points in Section 1. 

We can define zig-zag paths around white and black squares as in [8, 
9, 14, 15]. Consider the sequence of white squares in the Aztec diamond 
with opposite corners Q T k = (— r + k, n + 1 — k — r), k = 0, . . . , n + 1, where r, 
1 < r < n, is fixed. The zig-zag steps from Q r k to Q r k+ i go either one unit step 
to the right and then one unit step down (ES-step) or the other way around 
(SE-step), and in such a way that it does not intersect a domino. Similarly, 
we can define zig-zag paths around black squares between (— r, n — r) and 
(n — r, —r). It is not hard to see, compare the proof of Lemma 2.2 in [15], 
that the zig-zag paths around white squares have a red dot exactly when 
we have an ES-step around it, and that zig-zag paths around black squares 
have SE-steps around squares with green dots. In this way each zig-zag path 
corresponds to a unique particle configuration, and we can view the point 
process as the totality of all zig-zag particles. 

It was proved in [14] and [15] that the zig-zag particles along a single line 
define a point process described by the Krawtchouk ensemble. The possible 
positions of the red particles can be taken to be {0,1,..., n} and the rth 
zig-zag path has r red particles at h\ , . . . , h r (no ordering) . The probability 
of having particles at exactly these points is [15], Theorem 2.2, 



where q = a 2 (l + a 2 ) -1 , p = 1 — q and A r (h) is the Vandermonde deter- 
minant. This Krawtchouk ensemble has determinantal correlation functions 
given by the Krawtchouk kernel, 



(2.6) 



F[h] = -^—A r (h) 2 Y[(h) 1 hj P n ~ Hj 



(2.7) K r ^ q (x,y) = J2Pk(x;q,n)p k (y;,q,n) [ jffp' 




nl/2 



where 
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are multiples of the ordinary Krawtchouk polynomials; 7 is a circle centered 
at the origin with radius < min(l/p, 1/q). Pk(x; q, n) = (^,) 1 ^ 2 (q/p) k ^ 2 Kk(x; q, n), 
where K\. are the standard Krawtchouk polynomials, see, for example, [19]. 

Next, we want to show that the fact that we have a determinantal point 
process when we restrict to a single line can be extended to show that the 
whole point process is determinantal. Using the Lindstrom-Gessel-Viennot 
method, the combinatorial version of the Karlin-McGregor theorem (see, 
e.g., [26]) we can write the probability for a certain particle configuration x 
as a product of determinants 

j 2n-l 

(2.9) VN,n{x) = — II det (^,r+l(4',4 +1 ))jS=l' 

where the transition function T . jr+ i(x, y) gives the weight of all paths going 
from x on the vertical line r to y on the vertical line r + 1. The right- 
hand side of (2.9) is a symmetric function of x\,...,x r N for any r, so we 
need not consider the variables on a single vertical line as ordered. Prom 
the definitions of the directed graph Q' and its weights, we see that we can 
describe the nonintersecting paths using two types of transitions between 
different vertical lines. From vertical line 2i to 2% + 1, we can go from (2i,j) 
to (2i + l,j), with weight 1, or from (2i, j) to (2i + 1, j + 1), with weight a. 
From line 2i + 1 to 2i + 2, we have to take a step to the right, with weight 
1, and then a certain number > of down steps, with weight a. This leads 
to the following transition functions: 

]fy-X = l, 

(2.10) 4>2i,2i+i(x,y) = a(y - x) = { 1, ify-x = 0, 



and 

(2.11) (f>2 i+ i,2i+2(x,y) = (3{y - x) 



otherwise, 

a -(y~ x ), ify-x<0, 
0, otherwise. 



Set f2i(z) = az + 1 and f2i+i{z) = (1 — a/z) 1 so that a(n) = f2i(n) and 
(3{n) = f2i+i(n), where f r (n) is the nth Fourier coefficient of f r . We assume 
that < a < 1. The case a = 1 will be handled by taking a limit a — > 1— . Note 
that f r has a Wiener-Hopf factorization f r = f^~f^, where /^(z) = az + 1, 
f 2 ~i = fti+i = 1 an d f2i+i( z ) = (1 — a/z)' 1 . It follows from [16], Theorem 17, 
that the point process on V n = {1, . . . , 2n — 1} X Z defined by (2.9) has 
determinantal correlation functions. The probability of finding particles at 
z j = ( r j' %')' 1 — j — m > w here < r j < 2n and yj G Z, is 

(2.12) det(K Ntn (zj;z k ))™ k=1 , 
where the kernel is given by the formula 

(2.13) K NtTl (r,x;s,y) = K Nn (r,x;s,y) - (j> rs (x,y). 
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Here c/v,s = 4> r ,r+i * ■ ■ ■ * <j> s —i,s if s > r and r ,s = if r > s. Furthermore, 

TV 

(2.14) K Ntn (r,x;s,y)= ^ (f) r ,2n( x , ^ n )(^~ 1 )ij'hA x j' y), 

where A = (aij)f^ =l , aij = 4>o,2n{xi,%'j n )- This is a consequence of the fol- 
lowing formula, Proposition 2.1 in [16]. Assume that g : V n — > R is a bounded 
function. Then, 



(2.15) E 



7V,n 



2n-l JV 



det(I + gK Nt n) L 2(y n 



n n( i +^4')) 

r=l j=l 

where is counting measure on V n . Here Ejv,n is expectation with respect 
to the probability measure (2.9). If g(r, x) = 0, for x < —n + [r/2], then the 
left-hand side of (2.15) depends only on the part of the particle configuration 
that corresponds to the Aztec diamond. We know that this is independent 
of N, compare (2.3), and we can replace the kernel Kn,h d Y its limit as 
N — > oo. This limit is given by, Proposition 1.8 in [16], 

(2.16) K n (r,x;s,y) = K n (r,x;s,y) - <j> ryS (x,y), 
where 

— ' f dz f div vj z 

(2.17) K n (r,x;s,y) = / — / -G n%rs (z,w) 



(2ni) 2 J\z\=n z J\w\=r 2 w z x ' ' ' z — u?' 
a < r\ < 1/a, < r2 < r\, 

(2.18) Gn , rt8 (z,w) - nL -i /t+(1/z)n -- 1/r(1/u;) 
and 

(2.19) <M*,y) = ^ J\ l ^ e G n , r , s (e w ,e w )d9 

for r < s and <j) riS {x,y) = otherwise. 

Inserting the formulas for into (2.18) we see that, with £±,£2 G {0, 1}, 

(2 20) G (z in) - (l-«H-* +£2 (l + a/< 

(2.20) G n> _ £l , 2s _ £2 (z,U,) - {1 _ az)n -r + e l{1 + a/z) r ■ 

Using Cauchy's theorem, we can deform the integration contours to 7 given 
by w(t) = aie lt , —ir < t < ir, and T given by z(t) = a 2 + it, t S R, where 
< a.\ < a 2 < 1/a. We can then let a — > 1— to get the case a = 1 using the 
continuity of all expressions involved. Then, 

K n (2r,x;,2s,y) 

(2.21) 1 f dz f dw wy (I- w) n ~ s (1 + 1/ w) s z 



{2ni) 2 J r z Jj w z x (l - z) n - r (l + l/z) r z-w' 
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It follows that the probability of finding particles at positions Zj = (rj,yj), 
where yj > — n + [rj/2], 1 < j < m, is given by (2.12) with K n instead 
of KN n . Consequently, for all our computations we will use K n given by 
(2.16). We will call this kernel the extended Krawtchouk kernel. For future 
use we record the following consequence of (2.14). Let r\ < ■ ■ ■ < r m and 
let £±,...,£ m be given real numbers such that i k > —n + [rfc/2]. Then, for 
N>n, 



N,n 



(2.22) 



max xJ < £j , 1 < j < Tn 



Kk<N 



: E N . 



2n-l N 

\[(l+g{r,x])) 

r=l j=l 



det(I + gK n ) L 2 (Am)fl) , 



where g(r,x) = — l^ fe)00 )(x) if r = 1 < k < m, and g(r,x) = otherwise. 
Since x\ = maxi<fc<v x r k when the variables are ordered, we can combine 
(2.22) with (2.5) to study the boundary of the NPR. To use it to prove 
Theorem 1.1 we need some asymptotic results for the extended Krawtchouk 
kernel. These will be proved in the next section. 

To justify calling (2.16) the extended Krawtchouk kernel we should show 
how it is related to the Krawtchouk kernel defined by (2.7). Let be the 
normalized Krawtchouk polynomials given by (2.8) and set = if k < 
or k > n. Define, <q <1, p = l — q, 



(2.23) 



L n , q (r,x;s,y) = 



n 

k + s 



-1/2 



n 

k + r 



1/2 



x p k+r (x; q,n)p k+s (y; q,n)[w q (x)w q (y)} 1/2 



E 

fc=0 



n 
k + s 



-1/2 



n 
k + r 



1/2 



if r < s, and 
(2 24) L n, q {r,x;s,y) 

if r > s, where w q {x) = Q)q x p y ~ x . Note that L ntq (r,x;r,y) = K r ^ q (x,y), 
where K r ^ n ^ q is the Krawtchouk kernel given by (2.7). If q = a 2 (l + a 2 ) -1 , 
then 



xp k+r (x;q,n)p k+s (y;q,n)[w q (x)w q (y)} 1/2 



(2.25) 



K n (2(n -r) + l,x - r + l;2(n - s) + l,y - s + 1) 

-l-i 1/2 



(-l) r s L niq (r,x;s,y), 



where K n is the extended Krawtchouk kernel defined by (2.16)-(2.20). The 
prefactor in (2.25) cancels in all determinants and is not important. Hence, 
we can just as well think of L n ^ q as an extended Krawtchouk kernel. We will 
prove (2.25) in Section 5. 
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Let us briefly comment on some asymptotic results related to the behavior 
of the domino tiling of the Aztec diamond close to the point of tangency of 
the Arctic circle, that is, we consider the positions of the zig-zag particles 
on the rth line where r stays fixed as n — > oo. The corresponding limit for 
the normalized Krawtchouk polynomials is [19], 

(2.26) lim p k (qn + ^V2^;q,n) = {-l) k h k (07r 1/4 , 

n — >oo 

where h k (x) = 2~ k / 2 {k\)~ 1 / 2 'K~ l l 4: H k (x) are the normalized Hermite poly- 
nomials. Define an extended Hermite kernel by 

(2.27) K n , I (r,Z;s,r ] )= £ WO W^" (fW)/2 

j=— oo V V •>)' 

if r < s, 



(2.28) Kn,i(r,t;s, V ) 

if r > s. Here h k = if k < 0. Note that (2.27) is the ordinary Hermite kernel 
of GUE(r) if r = s. Taking the limit (2.26) formally in (2.23) and (2.24) gives 

lim (-n) s ~ r \/2npqL n , q (r, qn + ^V^npq; s, qn + r)\j2npq) 
= Kn,i{r,£; s,rj). 

The kernel i^H,i also occurs in the following problem. Consider a random 
matrix A from GUE(re), and let Ak be the h x k upper left corner of A, 
1 < k < n. The kernel i^H,i then describes the correlations between the eigen- 
values Xj(k, n) of Ak, 1 < j < k, 1 < k < n, which form a determinantal pro- 
cess. This follows with some work from [3] and [16]. For the Aztec diamond 
this leads to the result that the zig-zag particles along the lines 1, . . . , r, r 
fixed, n — > oo, behave like these eigenvalues, see [21]. Note that there is also 
a different extended Hermite kernel which describes the correlations of the 
eigenvalues of Dyson's Brownian motion at different times, namely, 

(2.30) K U:U (r,x;a,y)= ]T e k ^h n+k (x)h n+k (y) e -( x2+ y 2 V 2 

k=—oo 

if t > a, and 

oo 

(2.31) K U)n (r,x-a,y) = -Y / e k{T - C7) hn + k(x)h n+k (y)e-^ + y^ 2 

k=0 

if r < a. Again, (2.30) reduces to the ordinary Hermite kernel for GUE(ra) 
if r = a. 
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3. Asymptotics. The rescaled variables £, £',r and t' are defined by 

2r = 6 n (r) = n(l + 1/V2) + 2" 1/6 rn 2/3 = n(p + 1), 
2s = 6 n (r / ) = n(p / + l), 

x = n(p + Vp 2 -l)/2 + 2~ 5 / 6 ,en 1/3 = na(p) + 2" 5 / 6 1 en 1 / 3 , 

y = n (p' + VT^T )/2 + 2" 5 / 6 £V /3 . 

If we write a n (r) = - 2" 5 / 6 r 2 n 1 / 3 , then n(p + V p 2 - l)/2 = a n ,(r) + 

• • • . Set 

J^(2r,x;2 S ,y) 

(3 ' 1} = (V2- l)^ +2 ( S -'-) e ^ V -( 1 / 3 ^ 3 +( 1 /3)r' 3 ^ n(2r)X;2S)y)) 

where X n is given by (2.16). We can use _fC* as our correlation kernel instead 
of K n , since all determinants are unchanged. Set c n = 2 -5 / 6 n 1 / 3 and define 
the rescaled kernel 

(3.2) L* (r, £; r, £') = c n K*(b n (T) , a n (r) + c n £; & n (r'), a n (r') + c n £')- 

The basic asymptotic result is contained in the next lemma which will be 
proved below. 

Lemma 3.1. (a) Uniformly for ,r,r' in a compact set, 
(3-3) Jim L;(r^;r / ,0 = ^(e,T;e , ,r'), 

where A is the Airy kernel defined by (1.2). 

(b) Fix M > 0. There are positive constants d\,d2,N , which depend only 
on M such that if r > r', 

(3.4) iL^r^^'^Ol^^ie-" 2 ^ 
and if t < t' , 

(3-5) |^(r,ey,Ol < -^e-*^-^ 

V T — T 

/or a// n > AT, r, r' G [-M, M] and £, f ' E [-M, 00) . 

The convergence of Fredholm determinants will be proved using their 
Fredholm expansions. 

Lemma 3.2. Let t\ < • • • < r m and 71, . . . ,j m £ R be given. Set A m = 
{ri,...,r m } and let v be the counting measure on A m x Z. Also, write 
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g n {Tj,x) = -t(a n {T J )+- ii c n ,oo){x), xeZ. Then 

i™E^[ det{K n {z uZ] )) k l]=l \{g n {z 3 )d k u(z) 

(3-6) 

= Eu / A , det(A(^,^))t = iII/^) d V^), 

where [i = fi\ <&ii2, l^i "is a counting measure on A m , /12 is a Lebesgue measure 
on R and f(rj,x) = — l( 7j)0 o)( x ); 1 < J < m - 

Proof. The proof has many similarities with the proof of Lemma 3.1 
in [13], where more details can be found. The sum in the left-hand side of 
(3.6) can be written 



f^kl J(A m xZ)k C n V V Cn 

(3.7) T{z ^ X ( Zj )-a n (r( Zj )) Xy 

c n / / i,j=l 

X ]]/(rfe), ^^ )^), 

j=1 V ^ J 

where r(zj) = tj, x(zj) = Xj if Zj = (Tj,xj). Find M > so that Tj G [— M, M], 
7j > — M, j = 1, ... ,m. Then, by (3.4), (3.5) and the Hadamard inequality, 
the determinant in (3.7) is, ds,d± > 0, 

< 4 n ( e (-^4 (' w -^ (r(at)) + ^z^Mj ) ) 1/2 



i=l \j=l 

<4e feM ^A; fc / 2 nexp(-d4 



x(zi) - a n (r(zi)) 

for those Zj that contribute to the integral in (3.7). Using this estimate, we 
see that, given e > 0, the part of the sum in (3.7) where k > N is 

\ h 

e -d4,(x-a n {Tj))/cn I < £ 



00 -, / -, m 00 

<E^WfE E ■ 

fc=-/V ' \ n j=lx=a n (Tj)—Mc n 



•3 1 

if N is chosen large enough. Similarly, using the estimates (3.4), (3.5) and 
the Hadamard inequality again, we can restrict the integration in (3.7) to 
A m x ((— oo,Lc n ] n Z), with an error < e by choosing L large enough. We 
can then use (3.3) to see that what remains converges to 

AT . k 
( 3 - 8 ) Eu[ , n k ^m^ j )% =1 J{f{z j )d k Kz). 

k=0 • / (AmX(-oo,£])fc i= i 
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The extended Airy kernel satisfies estimates like (3.4) and (3.5), as can be 
proved using standard estimates of the Airy function. The same type of 
argument as above then shows that (3.8) approximates the right-hand side 
of (3.6), with an error < e provided N and L are chosen large enough. This 
completes the proof. □ 

We can now give the proof of Theorem 1.1. 

Proof of Theorem 1.1. From (2.22) it follows that, with 2rj = b n (rj), 



k < a n (Tj) + 7j-c n , 1 < j < m 



2r 

X 

l<k<N 

= det(I + gK n ) L2{Am xZ)I/) 



(3.9) 

"■ \ r * 
= Ep[ detiKni^z^l^Ylg^d^z), 

with g as in Lemma 3.2. Now the right-hand side of (3.9) converges to 

OO -1 „ k 

(3.io) zjsL^^^w™*"* 



k=0 "K^m^r j = l 

(A m xR)i 



det(I - fV*AfV\ H 



where f(rj,x) = 1( 7 )0 o)( x )- Combining (3.9), (3.10) and (1.3), we see that, 

N>n, 



, , lim F N n 

[6.LL) n— »oo 



max xP < n/V2 + 2~ 5/6 h i - T?)n 1/3 , 1< j < r 
i<k<N k ~ 1 v u 3 > ■> —j - 

P[^)<7j,l<j<H- 

To conclude the proof we use (2.5). Note that in (2.5) the variables are 

It ■ 1v It / — 

ordered so that x 1 3 = maxi<fc<7v x k J . Write x t 3 = n/v2 + c n Cj- Since 2rj = 
n(l + l/\/2 ) + 2" 1 / 6 r i n 2 / 3 , we get 

X n {2-^n^(r - 2^Q)) - + 1/2 

Combining (3.11) and (3.12), we obtain (1.4) and the theorem is proved. 
Note that the variation in the argument of X n due to Q is negligible. To 
be more precise, we could use the uniformity in t±, . . . ,r m in our convegence 
estimates. □ 

We still have to prove our basic asymptotic lemma. 

Proof of Lemma 3.1. Set 

F r , x {z) = z*- r (l-z) n - r (l + zy. 
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Then, by (2.21), 

^ ^ n( , , ,y) [2iri) 2 Iy z i-y w F T)X {z) z — w 

and, furthermore, 



(3.14) fe , l( ,,,) = -L/.»— (i_£)' 



l r /l — z\ r ~ s dz 

<p2r,2s(X,y) = - 

if r < s. Set 

/(*) = («(P) " ^) log z + ^ log(l - z) + log(l + *) 
so that 

(3.15) F r , x (z) = e 2 ' 5/ ^ nl/3lo ^ +n ^ z l 

A straightforward computation shows that f'(z) has a double zero z c = 
z c (p) = p(l + Vl — p 2 ) _1 . It is this double zero condition that specifies a(p), 
and corresponds to the arctic circle. The geometrical considerations and 
the choice of the place where we want to show convergence to the Airy 
process leads to p = l/y/2 + 2 _1//6 rn _1 / 3 . (The exact number 2" 1 / 6 comes 
out of the computations below.) We get z c = y/2 — 1 + r/dn 1//3 + • • • , where 
d = 2~ 5 / 6 (l + \/2). Also, we write z' c = z c (p'). 
As our paths of integration we will take 

z(t) = z c + j-^, t€R, 

w (t>) = U - -^-\jt/d>nV^ |^| < 7rd / n l/3 ) 

where d' = d(y/2 - 1) = 2~ 5 / 6 , and r? > is such that r - r' + 2i] > 0. 
Claim 3.3. 5ei 

Then with the choices ofr,s,x,y as above: 
(i) 

(3.17) lim A n ( V ) = e ^'-^+^/3-r' 3 /3 

n — >oo 

uniformly for ,t,t' ,rj in a compact set, and 

(ii) for any > —M, |r|, |r'| < M, < 77 < 3M, £/iere is a constant, 
which only depends on M , such that 

(3.18) |A n (ri)| < Ce ?V -^-"^'). 
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Proof. Write 

w c = z c + rj/dn 1/3 = \/2 - 1 + (r + rj)/dn 1/3 + ■■■, 

w ' c = z ' c - ri/dn 1 ^ 3 = y/2 - 1 + (r - rj) /dn 1 / 3 + ■■■. 

The higher-order contributions can be included in r, t' by slightly changing 
their values. Then, 

(V2-l) 2r - x F r!X (w c ) 

= w x c - r (i-w c ) n - r (i + Wc y 

= (1 + 2 5 / 6 (r + ^n- 1 / 3 )*-^/ 2 ^ - 2 1 / 3 (r + t 7 )^- 1 /3)— 

x (1 + (y/2- l)2 1 / 3 (r + rj)N-V 3 ) r . 

We can now insert x = n/y/2 + 2~ 5 / 6 (£ — T 2 )n l l 3 + • • • and r = n(l/2 + 
1 /2 \/2 ) + 2~ 7 / 6 rn 2 / 3 into this expression. Also, we can write an exactly anal- 
ogous expression for (y/2— l) 2s ~ y F Sjy (w' c ). A straightforward, but somewhat 
lengthy, computation now gives (3.17) and (3.18). □ 

Claim 3.4. Let r,s,x,y be as above and set 

n n = {(t, t')£Rx [-Trd'n^M'n 1 / 3 ]; \t\ > n 1 / 3 " 5 , \t'\ > n 1 / 3 " 5 } 
for some fixed 5 > 0. Then, uniformly for £, r, r' in a compact set, 

(3 - 19) ^y nB «,(to(z(t)-u,(f)) F r ,^(t)) )(V2 - 1} 

The claim will follow from the estimates used to prove Lemma 3.1(b) 
below. Let us accept it for the moment. Expanding f(z) around z c , we get 

(3.20) F r>x (z(t)) = e-^-^-WW-t+^+^^Fr^Zc), 

where rn\t) — > as n — > oo uniformly for \t\ < n 1 / 3 "* 5 and £, r in a compact 
set. Similarly, 

(3.21) F s Jw(t')) = e^'C'+^-^^+^^rfw^^^/^ 

where rffl (t) — > as n — ► oo uniformly for < n 1 / 3- ' 5 and r' in a compact 
set. If we insert this into (3.13) and use Claims 3.3 and 3.4, we obtain, after 
changing t to — t, 

lim (y/2- l) x ~ v+2{s - r h-^ & n 1 ' 3 K n (2r,x-,2s,y) 

n — >oo 

(l/3)(r 3 -r' 3 )-^+£V 



( 3 - 22 ) - ^ e i£(t+ir,)+iti'(t'+iv)+i((t+iv) 3 +(t'+iv) 3 )/3 

dt / dt'- 



t' -r + i(t + t' + 2rj) 
eW 3 ^-T' 3 )-^+^'A(T,t-y,t'), 
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by Proposition 2.3 in [16]. In order to prove Lemma 3.1(b) and also Claim 
3.4, we need some estimates. 
Write 

(3.23) \FrAz(t))\ = \z(t)\ 2 ~ 5/6 Z niri \F r , X0 (z(t))\, 

where x = na(p) = n/y/2 — 2 _5 / 6 r 2 n 1 / 3 H . Now, 

\F r , X0 (X + iu)\ 2 = (A 2 + u 2 T»~ r ((l - A) 2 + u 2 ) n - r ((l + A) 2 + u 2 ) r 
and, hence, 

±\ og \F r , X0 {\ + iu)\ 2 = 2u? {u) 



du q{u) 
where 

p(u) = (x - r) [(1 - A) 2 + u 2 ] [(1 + A) 2 + u 2 } 

+ (n - r) [A 2 + u 2 ] [(1 + A) 2 + u 2 ] + r [A 2 + u 2 } [(1 - A) 2 + u 2 ] 

and 

q (u) = (A 2 + u 2 )({l + A) 2 + u 2 )((l - A) 2 + u 2 ). 

We can now insert (x — r)/n = —1/2 + 1/2 \/2 — 2~ 7 / 6 rn~ 1 / 3 + • • • ,r/n = 
1/2 + 1/2 v/2 + 2- 7 / 6 rn- 1 / 3 and A = a/2 - 1 + 2" 5 / 6 (a/2 - l)(r + ^n -1 / 3 
into After some computation, we get 

-p(u) = 2 4 / 3 ( v / 2 - 1) V~ 1/3 + (V2 - l)u 2 + (l/l/v^u 4 , 
n 

up to negligible contributions. Assume first that < u < 1; the case — 1 < 
u < is analogous. Then there is a numerical constant cq such that 

^log\F r>X0 (X + iu)\ > 2c r/re 2/3 u, 



and, consequently, 



\F T ^(X + iu)\>\F r , X0 {\)\e c ^ u \ 



for |tt| < 1. (Below cq denotes a positive numerical constant the value of 
which may change.) This gives, \t\ < dn 1//3 , 

(3-24) \Fr,x (z(t))\ > \F rtXo (z(0)) | e c °< 

If -u > 1, we use p(u) > nu 4 and g(u) < cou 6 and, hence, 

^log|-F r)a . (A + m)| >c n/u. 



20 K. JOHANSSON 

Combining this with (3.24), we obtain 



(3.25) 



\Fr,x o (z(t))\>\F rtX0 (z(0))\e Co " n 



2/3 



dn 1 / 3 



for \t\ > (in 1 / 3 . Looking back at (3.23), we see that we must also consider 



(3.26) 



zit) 



z(0) 



2-5/6 ?n l/3 



3 2- 1 /6(2 ?n -l/3 + . 



If |£| < M, the factor (3.26) can be absorbed into the other factors in the 
estimates above by slightly changing the constant c$. If £ > M > 0, then the 
expression in (3.26) is > 1 and we use just this trivial estimate. Combining 
the estimates, we obtain 



(3.27) 
where 



\F r Mt))\ > \F r>x (z(0))\g n (t), 



9n(t) 



exp(coryi 2 ), 

exp(c r?n 2 / 3 ) \t/dn 1/3 \ c ° n , 



if \t\ <dn 1 / 3 , 
if \t\ > dn 1 / 3 . 



We also have to estimate the w-integral. Consider 



\F S:y (w)\ = \w\y~ s \i 



w 



\l + w\ 



set w = ae td and define g{6) = (re — s) log |1 — ae ld \ 2 + slog |1 + ae lS \ 2 . Then, 

.(1 - 2/3) (1 + a 2 ) + 2acos<9 



(3.28) 



g'(0) = 2ansin6»- 



(l + a 2 ) 2 -4a 2 cos 2 (? ' 

where 2(3 = 2s /n = 1 + 1/^2 + 2~ 1 / 6 r / n~ 1/3 . The numerator in (3.28) is 
< for all if a < z' c = y/2 - 1 + 2" 5 / 6 (y2 - l^n" 1 / 3 . Hence, this will be 
satisfied if we choose 



dn 1 / 3 



Jt'/d'n 1 / 3 



with n > as above. A computation gives g(9) — g(0) < —co9 2 n 2 / 3 ri and, 
thus, 



(3.29) 



F. 



dn 1 /' 3 



-C t l2 Tj 



for \t'\ < ird'n 1 / 3 . Combining the estimates (3.27) and ( 3.29), we get 
\K n (2r,x;2s,y)\ 



< 



C IF, 



M-ri/(dnV*))\ ( [°° 

f v J_«/M«l/3^l / .>W 



n^lF^^ + nAdnVS))! V 



Trd'n 1 / 3 
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where we have used (3.18) in the last inequality. Let rj = 2t]q. Then, £'(t' — 
v) ~ £( r + V) ^ ~Vo(£, + f° r £' — 0) if we choose ryo so that t' — t]q < 0, 
r + % > 0. This proves (3.4) for the K part, and we also use the same 
estimates to show Claim 3.4. Notice that r 1 — r — 2t] < if t]q is chosen large 
enough. 

It remains to consider the 4> r ,s(x,y) part of K n in (2.16). Let r,x,s,y be 
as above with r' > r. We want to show that, for £,£',t, t' in a compact set, 
we have 

Cn (V2-l) x ^ +2( ^ r V 2r ,2 S (x,y) 

(3.30) _^ 1 p -(€ , -g+r 2 -r' 2 ) 2 /4(r'-r) 

^JAtt(t' - t) 
uniformly as n —> oo, and that, given M > 0, 

(3.31) |c n (V2 - l)^ y+2(s - r) ^-^V 2r , 2s (x, y )| < -^E^W-^-K') 

VT — T 

for |r|, |t'| < M and re sufficiently large (depending on r' — r). 
To prove this we use the integral representation 

(3.32) ^rt.-^^r-^Qli)", 

where < a < 1. Set z = ae ie and /(0) = logfll + ae ie \ 2 \l - ae i9 \~ 2 ). Then 

1 4- a 2 

/'(0) = -4a sin0- 5 ±f 5 -, 

v ; (l + a 2 + 2acos<9)(l + a 2 -2acos<9) 

so we have a maximum when 9 = 0, and 

(3.33) - /(0) < -c # 2 

for |0| < 7T. We can use this estimate to localize the integral to a small 
neighborhood of 9 = and a standard argument then proves (3.30). Prom 
(3.33), we obtain 

/ 1 _L ry\ s-r 



l-a) ?T, 1 / 3 v /47rco(r / - r) 
Hence, 

C n ( V2 - lf-V+«-r) e ir-i V | ^ (Xj y ) | 
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Choose a = (\/2 — 1)(1 + 2~ 1 / 6 n~ 1 / 3 (T + r')). Inserting the expressions for 
x,y,r,s into (3.34), a computation shows that the right-hand side of (3.34) 
is 

< c ° e -(T'-T)(g+g')/2+cm-y 3 (|g|+|g'i) 

yV — r 

for |r|, |r'| < M, where ci depends on M. Hence, if n is large enough, this 
expression is 

since r' — r > 0. This completes the proof of Lemma 3.1(b) and, hence, of 
the whole lemma. □ 

4. The center of the Aztec diamond. It has been conjectured by Propp, 
see [5, 6, 12], that in the limit n — > oo, a random domino tiling of the Aztec 
diamond "looks like" a random domino tiling of the whole plane under the 
Burton-Pemantle measure (the tiling measure with maximal entropy) [4] . In 
this section we will discuss and outline a proof of a version of this conjecture. 
The version given below is in terms of particle configurations and not directly 
in terms of dimers. A domino tiling of the Aztec diamond or the whole plane 
can equivalently be described as a dimer covering (perfect matching) of a 
certain graph. Let V be the set of vertices (1/2 + j, 1/2 + k), j,k £ Z, in 
CS-I and E the set of edges between nearest neighbors. A domino tiling 
of the whole plane is the same as a dimer covering of (V, E). It will be 
convenient to write the coordinates as complex numbers. Color the point 
1/2 + j + i(k + 1/2) black if j + k is even and white otherwise. Let w he, a, 
white vertex in V and give an edge between w and w + z, where z = ±1, dbz, 
the weight z. Let D C E be a finite subset of edges in E and assume that 
the edges in D cover the black vertices b\ , . . . , b m and the white vertices 
w\, . . . ,w n . It is proved in [18], using techniques going back to Kasteleyn [17], 
that under the Burton-Pemantle measure fi, the probability of the event Ud 
that the edges in D belong to the dimers of a dimer covering of (V, E) equals 

(4.1) v(U D ) = a D det(P(b j -w k ))™ k=1 , 
where is the product of the weights of the edges in D and 

1 r w ^(xQ-yfy) 

(4.2) p( x + iy) = — / -d6d(j). 

v ' \ y> 47r 2 y_ ff< /_ w 2isin0 + 2sin0 Y 

Let A n be the Aztec diamond region as before, set V n = A n n V and let E n 
be the edges between nearest neighbors in E n . Then dimer coverings of the 
graph (V n ,E n ) are in one-to -one correspondence with domino tilings of A n . 
Assume that n is odd so that the square with center 1/2 + j + i{k + 1/2) 
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is black if and only if j + k is even. The case of even n is analogous. The 
red and green particles defined in Section 1 lie in V n . It follows from their 
definition that there is a red particle at v if and only if v is white and there 
is a dimer from v — 1 to v or from v — % to v. There is a green particle at v 
if and only if v is black and there is a dimer from v to v + 1 or from v to 
v + i. The discussion in Section 2 also shows that these particles uniquely 
determine the domino tiling of A n . We can similarly associate a red/green 
particle process with domino tilings of the whole plane by using the relation 
to dimers just described. From (4.1) we can then compute the probability 
of seeing red/green particles at specified points. Note that all particles on 
the line y\ = —x\ + 21 + 1 have to be green, whereas all the particles on 
the line y\ = —x\ + 2£ have to be red. Consider, for simplicity, only green 
particles. We wish to compute the probability of seeing green particles at 
Vj = Uj + 1/2 + i{—Uj + 2£j + 1/2), 1 < j < m. Then all Vj are black and 
there has to be a dimer between Vj and Vj + 1 or between Vj and Vj + i. 
Hence, this probability is given by, compare Theorem 2.10 in [15], 

zi---z m deb(P(vj - (vj + zfc)))™fc=i 

Zj=l or i 

m 

«i ■ • -*m s S n ( £J ) II P ( v i ~ Mi) ~ Mi)) 

Zj = l or i creS m J=l 

m 

(4.3) sgn(a) ^ J] z *(j) P ( v i ~ v *ti) ~ Mi)) 
J2 sgn(a) f[l J2 P ( V J ~ Mi) ~ Mi) ) J 

a-eS m j=l \z ff (j)=l or i / 

= det(P(vj -v k -l) + iP(vj -v k - i))™ k=v 

Set, v G C, 

(4.4) R(v) = P(v-l) + iP(v-i). 

Consider now the probability of seeing green particles at v i , . . . , v m in the 
Aztec diamond A n , with n odd and sufficiently large to include these points. 
This probability can be expressed in terms of the extended Krawtchouk 
kernel. A green point corresponds to a point Q + (1/2, 0) in CS-II where Q is 
a last particle position. In CS-I this corresponds to a position, compare (2.5), 
(2r-n + l/2-x 2r ,x 2r -l/2), which gives 2r = 2£ + n + l, x 2r = 2£-u + l. 
Hence, in terms of the extended Krawtchouk kernel, (2.16), the probability 
of finding green particles at V\ , . . . , v m is 

(4.5) det(K n (2£j + n + l,2£- Uj + 1; 2£ k + n+l,2£ k - u k + l))™ k=v 
We can interpret Propp's conjecture as saying that 

/ lim det{K n (2£i + n+ 1,2£- Uj + l;2£ k + n + 1,24- u k + l))™ k _ x 

(4 6) n^oo J ' ' 

= det(R(uj -u k + i(u k - uj + 2(£j - £ k ))))J k=1 . 
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The proof of (4.6) consists of two steps. First, we must compute the asymp- 
totics of the kernel K n for the appropriate values of the arguments corre- 
sponding to the center of the Aztec diamond. We will do this using the 
approach in [22], but we will not give the details, which are of a rather stan- 
dard saddle-point argument nature. We then have to show that the resulting 
limiting expression equals the right-hand side of (4.6). These computations 
will be presented below. The cases of just red particles or a combination of 
red and green particles can be treated completely analogously. 
It follows from (2.16), (2.21), (2.19) and the residue theorem that 

K n (2r,x;2s,y) 

(4.7) 1 f dz f dw wy (I -w) n ~ s (1 + 1/ w) s z 

~ (2ni) 2 Jr ~~z ~w z x (l - z) n ~ r (l + l/z) r z-w" 

with r : t — > ct2+it, t E M, 7: t — * aie lt , \t\ < tt, where 02 > «i > if r > s and 
a.2 < — ai, cti > if r < s. Write n = 2N — 1, Xj = 2£j — Uj + 1, 1 < j < m. 
Then, 

K2 N -i(2(£j + N), Xj ; 2(4 + N),x k ) 

(4.8) 1 f dz f dw w Xk {l-w) N (1 + 1/ w) N+tk z 

~ (2m) 2 Jr~ J-y ~w z x i(l- z) N ~ l -^(l + l/z) N+ ^ z-w' 

Write (1 - w) N (l + l/w) N = w~ N (l - w 2 ) N = exp(NF(w)), where f(w) = 
— logw + log(l — w 2 ). Then, f'(w) = if and only if w = +i. Assume that 
lj > 4; the other case is similar. We can move T so that it goes through +i, 
TqU — > it. We then pick up an extra contribution from the pole at z = w 
when w is on 7+ : t — > e , —tt/2 <t< ir/2. We obtain 

K 2N -i(2(£ j + N), Xj ;2(e k + N),x k ) 

(4.9) 2vrii 7+ w v ' v 11 

1 r dz f dw w x « (1 -w)"- 1 "^ (1 + 1/ w) N+£k z 
+ (2vri) 2 7r z J 1 V z x j(1- z) N - 1 -^(l + l/z) N+e J z-w' 
The second integral goes to zero as iV — > 00, compare [22], and we obtain 
lim K 2N -i(2(l 3 + N),2£j - Uj + 1; 2(4 + N),2£ k -u k + l) 

(4.10) =± r dw^_ u J 1-w 

27TZ J 1+ W \w(l + w) J 

It follows from (4.4) and (4.10) that in order to prove (4.6) it suffices to 
show that 

(AH) 2m h+ w \w(l + w)J 

1 ' ' =i u - v+2 ^- m ^[P(u-v-l+i(v-u + 2(£-m))) 

+ iP(u-v + i(v -u- 1 + 2(1- m)))]. 
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We will use the following formula. If x + y > 1 and x + y is odd, then 



(4.12) 



P(x + iy) 



2vri 



7+ 



f3! - B+ i)/2 (ti>-l) (3C+y - 1)/2 

{w + l){x+y+l)/2 w ■ 



Assume that we have proved (4.12). If we set x = u — v — 1, y = v — u + 2( 
m), then the right-hand side of (4.11) is 



.(_l)(*+rt-i)/2^+i J_ / ( 



1 

2vri 



IV 



2m 

(x-y+l)/2 



[{ x-l +j x+l w j w (x-y+l 



)/2 (w - i)(*+y-i)/2 dw 
(w + i)(s+y+i)/2 „, 



7+ 



W + iy 



n (x+y+X)/2 dw 



W 



which is exactly the left-hand side of (4.11). It remains to prove (4.12). 
Make the shift 9 9 - (f) in (4.2) to get 

e -i(x+y)tp 



(4.13) P(x + iy) 
Since 2is'm(9 — 6) 



1 

~ 4^ 
+ 2sin 

1 

~ 2^ 



tt 2isin( 



+ 2 sin c 



e ix9 d9. 



[e id + i - {e~ ie + i)e 2 ^]e"^, we obtain 



= i((:r+y-l)/2)20 



tt e ie + i 



1 + i)e 2i< ^ 

-(x+y-l)/2-l 



d9 



dz 



d9, 



_2m' J\z\=i 1 — az 
where a=(e- ie + i){e ie + If we make use of the integral, compare [22], 



1 

2ni 



-k-l 



z=l 



1 — az 



■ dz 



0. 



if k > 0, |a| < 1, 
if k < 0, |a| > 1, 
otherwise, 



we find 



P(x-Hy) = — 



„ix9 



o e* 61 + % V e ?y + i 



-i0 



+ i 



(x+y-l)/2 



de, 



which is easily seen to equal the right-hand side of (4.12). A similar com- 
putation should relate the Kasteleyn kernel for the finite Aztec diamond, 
computed in [11], to the extended Krawtchouk kernel K n . 

5. Extended Krawtchouk kernel and Krawtchouk polynomials. In this 
section we will prove (2.25), which expresses the extended Krawtchouk kernel 
in terms of Krawtchouk polynomials. We have the formula 

K n (2r - ei,x;2s - e 2 ,y) 



(5.1) 



dz 



dw w y (1 — aw) n s+S2 {l + a/w) s z 



(27ri) 2 7 7ri z J 7r2 w z x {l-az) n - r+e ^{l + a/z) r z-w' 



26 K. JOHANSSON 

where £i,£2 £ {0, 1}, a < r\ < 1/a and < r%,ri if 2r — e\ > 2s — £2, < r\ < 
r2 if 2r — £1 < 2s — £2- Here 7 r denotes a circle of radius r around the origin. 
This follows from (2.16)-(2.20) and the residue theorem. 
Consider first the case r < s. Then 

K n (2{n - r) + 1, x - r + 1; 2(n - s) + 1, y - s + 1) 
(5.2) _ (-l) s - r f dz f dw z n ~ x {aw - l) s {w + a) n ~ s+1 z 

(2vri) 2 i 7ri T7 7r2 ~w w n ~v (az-l) r (z + a) n - r+l z-w' 

Since < \z\, we have 

^ ((w + a)/{aw-l)) k _ (aw-l)(z + a) 
t'o ((z + a)/(az-l)) k ~ ( a 2 + l)(w-z) 

and, hence, the integral in the right-hand side of (5.2) can be written 
(-l) s - r (a 2 + 1) ™{ [ z n ~ x dz\ 



(2^)2 ^ V7 7ri (az - iy~ k (z + z J 



x 



(aw - ly-^iw + a) n+k ~ s+1 dw 

w n ~y w 



In the z-integral we make the change of variables z = (1 — oC)(C + a ) *j C = 
(l-az)(z + a)~ 1 . This maps 7 n to a circle surrounding the origin but with 
the opposite orientation. We can deform this circle to 7 ri using Cauchy's 
theorem. Hence, (5.3) equals 

1 g 1 f / (l-aO^XC + a)^ 

(5.4) 



(27ri)2^ o ( a 2 + l)«VA n C- fc - x C 



x 



(1 - auy- k - 1 (u> + a) n+fc " s+1 
. 2 ~ 



Here we have also used the fact that, by Cauchy's theorem, the ^-integral 
vanishes if k > r. We now set lo = w/a, C = &(1 + a 2 )~ 1 z and use p = (1 + 
a 2 ) -1 , ^ = a 2 (l + a 2 )~ 1 . A computation, where we also replace k by r — 1 — k, 
shows that (5.4) equals 



(5.5) 



"(1 + a 2 y- s a x ~y yV /" (1 ~ ^) n ~ x (l + ps)* dz 



(1 - w) s - r+k (pw + g )^-^- fc dw 

w n ~ y w 



By Cauchy's integral formula, this is the coefficient of w n v in 

r-l 



l + a 



-) S r a x - y Y(l - w) s " r+k (pw + q \n+r-s-kJ_ 



k=0 

(1 — qz) n ~ x (l +pz) x dz 



X 



7n 



Z* 1 Z 
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This can be written 



1 + 



\ x - y Y, {l-w) k {pw + q) n - k — 

k=s—r 



2m 

(l-qz) n - x (l+pz) x dz 

z 



y k+r- 



1 + 



-1 n 



1 



(5.6) 



a *-v ^2 YW-w) k {pw+ q r- k — 

(1 - qz) n ~ x : (1 + pz) x dz 
z 



7 n 



S 



1 + a 5 



a x ~ y 



k=s—r 



-1 



(pq)~ 



i 



27TZ 7 7r2 CJ fc+1 



I/=0 



x£(")(l- w) y {pq) y U y (pw + q) n - 



(1 — qz) n x (l+pz) x dz 

z 



yk+r—s 



By the binomial theorem, the y-sum equals 

(pqcu — pqwu + pw + q) n = {pquj + q + w(p — pquj)) n 



Y[ n v )^ + q) y {p-pqu J ) n - y w n -y. 



Hence, the coefficient of w n y in (5.6) equals, by (2.8), 



k=s—r 



1 + a 2 

1/2 



x {-l) k {pq) k ' 2 [I) Pk(y,q,n)(-l) k+r - s (pq^ k+r - s y 2 



n 

k + r — s 



n\ n 
y) \x 



1/2 



p k+r -s(x;q,n) 



s-l 

x E 



1/2 



-1/2 / n X 1/2 



fc=s— r 

which proves (2.25) when r < s. 



k + r — s 



p k+r - s (x;q,n)p k (y;q,n) 
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The case when r > s will be handled by reducing to the previous case. In 
the integral in 

K n (2(n - r) + 1, x - r + 1; 2(n - s) + 1, y - s + 1) 



■1)' 



(2vri) 5 



dz 



dw z r ' 



(aw — l) s (w + a) 



n—s+l 



; 7l . 2 w w n y (az — l) r (z + a) n r+1 z — u>' 

with o < r\ < 1/a, r<i > r%, we make the change of variables w — ► —1/w, 
z—>—l/z. If r3 = 1/n , ^4 = 1/^2 , we see that 

K n (2(n - r) + 1,3 - r + 1; 2(n - s) + l,y - s + 1) 



{2m) 2 
(1 - ai«) n - s+1 



dz 



Tr 3 



>4 



^ w n-(n-y) 



re-s+l)+l 



(1 - az) n - r+1 (z + a ) n -( n - r+1 ) +1 z-w 
= -(-l) v ~ x+r ~ s K(2(n - (n - r + 1)) + l,n - a; - (n- r + 1) + 1; 

2(n - (n - s + 1)) + l,n - y - (n - s + 1) + 1). 
By our previous computation this equals 

-l-i 1/2 



-(-1) 



y—x+i — s 



-1/2 



n 

k + s — r 



1/2 



(5.7) 



E 

k=r— s 

x p k+s _ r (n - x;q,n) 
x Pk(n-y,q,n) 
x [io g (n - x)w,(n - y)] 1/2 

Using the integral formula, (2.8), for the normalized Krawtchouk polynomi- 
als, it is not difficult to show that 



(5.8) 



Pk (n - x;q,n) = (-l) fc " V /2_ V n/2+ V-fc(x; <?, n). 



If we use this formula in (5.7) we obtain, after some simplification, 



1/2 



(-ir 



: E 

k=r— s 



-1/2 



n 

k + s — r 



1/2 



1/2 



-ir s E 

j=0 



x Pn-k(y,q,n)[w q (x)w q (y)] 1/2 

-1/2 / n x 1/2 



n 
j + s 



n 
j + r 



xp j+r (x;q,n) 

xPj+s(y,q,n)[w q (x)w q (y)} 1/2 , 
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where we have put j = n — k — s and extended the summation using pk = 
if k < or k > n. This is what we wanted to prove. 

Acknowledgments. I thank A. Borodin for a helpful discussion concern- 
ing extended kernels and E. Nordenstam for help with the pictures. 
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